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We summarize three recent efforts to constrain the first few moments of cosmic creation 
before and during the epoch of inflation. We consider two means to explain a slight dip 
in the power spectrum of the cosmic microwave background for multipoles in the range 
of I = 10 — 30 from both the Planck and WMAP data. We show that such a dip could 
be the result of resonant creation of a massive particle that couples to the inflat on field. 
For best-fit models, the epoch of resonant particle creation reenters the horizon at wave 
numbers of fc* ~ 0.00011 d= 0.0004 (h Mpc -1 ). The amplitude and location of these 
features correspond to the creation of a number of degenerate fermion species of mass 
15/A 3 / 2 m p i during inflation where A is the coupling constant between the inflaton 
field and the created fermion species. Alternatively, one can explain the existence of such 
a dip as due to a jump in the inflation generating potential. We show that such a jump 
can also resolve the excessively large dark flow predicted from the M-theory landscape. 
Finally, we summarize our efforts to quantify constraints on the cosmic dark flow from 
a new analysis of the Type la supernova distance-redshift relation. 
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1. Introduction 

Analysi£E of the power spectrum of fluctuations in the cosmic microwave back¬ 
ground (CMB) provides powerful constraints on the physics of the very early uni¬ 
verse. The most popular account for the origin of the primordial power spectrum is 
based upon quantum fluctuations generated during the inflationary epoch. 3 4 ] Sub¬ 
sequently, acoustic oscillations of the photon-baryon fluid distort this to produce 
the observed features in the angular power spectrum of temperature fluctuations in 
the CMB and the spatial power spectrum of matter density fluctuations. 

In this context, there now exists the highest resolution yet available in the de¬ 
termination of the power spectrum of the CMB from the Planck Satellite!®] In 
Mathews et alP we made note of a peculiar feature visible in the power spectrum 
near multipoles £ = 10 — 30. This is an interesting region in the CMB power spec¬ 
trum because it corresponds to angular scales that are not yet in causal contact, so 
that the observed power spectrum is close to the true primordial power spectrum. 

An illustration of the Planck observed power spectrum in this region is visible on 
Figure [l] from Mathews et al. 5 Although the error bars are large, there is a noticeable 
systematic deviation in the range £ = 10 — 30 below the best fit based upon the 
standard ACDM cosmology with a power-law primordial power spectrum. This 
same feature is visible in the CMB power spectrum from the Wilkinson Microwave 
Anisotropy Probe ( WMAP ),® and hence, is likely a true feature in the CMB power 
spectrum. This can be interpreted as an artifact of the cosmic variance,® or a phase 
transition in the inflation-generating potential. 7 9 However, here we adopt the the 
premise that this could be a real feature in the primordial power spectrum produced 
by new trans-Planckian physics occurring near the end of the inflation epoch. In 
particular, we show that this feature is well represented by the resonant creatiorG-21! 11 ! 
of Planck-scale particles that couple to the inflaton field as shown by the solid line 
in Figure [l] and we now describe in detail. 

This interpretation has the intriguing aspect that, if correct, an opportunity 
emerges to use the CMB as a probe of physics at the Planck scale ( m p i ~ 10 19 
GeV). Indeed, massive particles generically exist in Planck scale compactification 
schemes of string theory from the Kaluza-Klein states, winding modes, and the 
massive (excited) string modes. We have explorecP the possibility that the mass of 
some Planck-scale particles lies in the range of the vacuum energy associated with 
the last 9-10 e-folds of inflaton accessible to observation. 

2. Inflation Resonant Particle Production 

The details of the resonant particle creation paradigm have been explained in 
Refs. [5][T0j[lTj Here, we summarize essential features. In the basic inflationary pic¬ 
ture, a rapid early expansion of the universe is achieved through the vacuum energy 
from an inflaton field. In the minimal extension from the basic picture considered 
here, the inflaton is postulated to couple to massive particles whose mass is of order 
the inflaton field value. These particles are then resonantly produced as the field 
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Fig. 1. CMB power spectrum in the range of i = 0 — 50 from Mathews et al.E] Points with error 
bars are from the Planck Data Release® The dashed line shows the best standard ACDM fit to 
the Planck CMB power spectrum based upon a power-law primordial power spectrum. The solid 
line shows the best fit for a model with resonant particle creation during inflation. 


obtains a critical value during inflation. If even a small fraction of the inflaton field 
is affected in this way, it can produce an observable feature in the primordial power 
spectrum. In particular, there can be either and excessJSH 11 or a dip^in the in the 
power spectrum. Such a dip, will occur at an angular scale corresponding to when 
the epoch of resonant particle creation crossed the Hubble radius during inflation. 

We note that particle creation corresponding to an imaginary part of the effective 
action of quantum fields has been considered in.^In that case the same creation 
should occur at the present time. In that case compatibility with the diffuse 7 -ray 
background can be used to rule out the possibility of measurable effects from this 
type of trans-Plankian particle creation in the CMB anisotropy. However, the effect 
of interest here is a perturbation in the simple scalar field due to direct coupling 
to Plack-mass particles at energies for which the inflation potential is comparable 
to the particle mass and cannot occur at the present time. The present scenario, 
therefore is not constrained by the diffuse gamma-ray background. 

In the simplest slow roll approximation 3 “5 for the generation of density perturba¬ 
tions during inflation, the amplitude, <5#(fc), of a density fluctuation when it crosses 
the Hubble radius is just, 

, (1) 

07 Up 

where H is the expansion rate, and <fi is the rate of change of the inflaton field when 
the comoving wave number k crosses the Hubble radius during inflation. If resonant 
particle production affects the inflaton field, then the conjugate momentum in the 
field <fi is altered. This causes either an increase or a diminution in S H (k) (the 
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primordial power spectrum) for those wave numbers which exit the horizon during 
the resonant particle production epoch. 

For the application here, we adopt a positive Yukawa coupling of strength A be¬ 
tween the inflaton field (p and the field ip of N fermion species. This differs from 1 -^ 11 
who adopted a negative Yukawa coupling. With our choice, the total lagrangian 
density including the inflaton scalar field (p, the dirac fermion field, and the Yukawa 
coupling term is then written, 

Act = ^d^cpd^tp- V{(p) 

+ i/p^ip — mipip + NXcpipip . (2) 

For this Lagrangian, it is obvious that the fermions have an effective mass of 

M(cp) = m — NXcp . (3) 

This vanishes for a critical value of the inflaton held, = m/NX. Resonant fermion 
production will then occur in a narrow range of inflaton held amplitude around 

(p = (p*. 

As in Refs. mmm we label the epoch at which particles are created by an 
asterisk. So the cosmic scale factor is labeled a* at the time at which resonant 
particle production occurs. Considering a small interval around this epoch, one can 
treat H = H * as approximately constant (slow roll inflation). The number density 
n of particles can be taken as zero before t* and afterwards as n = ra*[a*/a(t)] 3 . 
The fermion vacuum expectation value can then be written, 

{ipip) = ?t*0(t — f*)exp [—3i?*(t — £*)] . (4) 

where 0 is a step function. 

Then following the previous derivations, 5 10 11 we have the following modihed 
equation of motion for the scalar held coupled to ip: 

'<P + 2,H/} = -V\(P) + NX{/pip) , (5) 

where V'{(p ) = dV/dcp. In the slow roll approximation, (</> ~ 0, H ~ Ff*, V'(cp) ~ 
V'(cp*) this differential equation after particle creation (£>£*) is then the same as 
iiUiULU but with a sign change for the coupling term, i.e. 

cp(t >£*) = </>* exp [—3 H(t - £*)] 

- t 1 ~ ex P h 3H (. f ~ **)]] 

NX 

+ ^- 5 -n*(£-£*)exp[-3iF*(£-£*)] . (6) 

Oil 5|s 

The physical interpretation here is that the rate of change of the scalar held rapidly 
increases due to the coupling to particles created at the resonance <p = (p*. 

Then, using Eq. (JTJ) for the huctuation as it exits the horizon, and constant 
H « H„ in the slow-roll condition along with 


din a = Hdt , 


(7) 
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one obtains the perturbation in the primordial power spectrum as it exits the hori¬ 
zon: 

S H = - [SH(a)] N x=o - . (8) 

1 + 0(a — a*)(lVAn*/|</>»|I7*)(a*/a) 3 In (a/a*) 

Here, it is clear that the power in the fluctuation of the inflaton field will diminish as 
the particles are resonantly created when the universe grows to some critical scale 
factor a*. 

Using k*/k = a*/a, the perturbation spectrum Eq. Q can be reducedMD] to a 
simple two-parameter function. 


S H (k) = 


_[£g(a)]ATA=0_ 

1 + 0(fc — fc*)A(fc*/fc) 3 In (fc/fc*) 


( 9 ) 


where the amplitude A and characteristic wave number &* (k/k t > 1) can be fit to 
the observed power spectra from the relation: 


fc* = — , (10) 

rus 

where ri ss is the comoving distance to the last scattering surface, taken here to be 
14 Gpc. 

The connection between resonant particle creation and the CMB temperature 
fluctuations is straightforward. As usual, temperature fluctuations are expanded in 
spherical harmonics, ST/T = a imYim(0, 4>) (2 < l < oo and — l < m < l). 

The anisotropies are then described by the angular power spectrum, Ci = (|a; m | 2 ), 
as a function of multipole number l. One then merely requires the conversion from 
perturbation spectrum Sn(k) to angular power spectrum C). This is easily accom¬ 
plished using the CAMB codef^ When converting to the angular power spectrum, 
the amplitude of the narrow particle creation feature in <5#(fc) is spread over many 
values of i. Hence, the particle creation feature looks like a broad dip in the power 
spectrum. 

We have made a multi-dimensional Markov Chain Monte-Carlo analysi£S 15 of 
the CMB using the Planck datcf 1 and the CosmoMC code ? 5 For simplicity and 
speed in the present study we only marginalized over parameters which do not alter 
the matter or CMB transfer functions. Hence, we only varied A and fc*, along with 
the six parameters, flbh 2 ,fl c h 2 ,9,T,n s ,A s . Here, fib is the baryon content, fl c is 
the cold dark matter content, 6 is the acoustic peak angular scale, r is the optical 
depth, n s is the power-law spectral index, and A s is the normalization. As usual, 
both n s and A s are normalized at k = 0.05 Mpc -1 . 

Figure [ 2 ] shows contours of likelihood for the resonant particle creation parame¬ 
ters, A and fc*. Adding this perturbation to the primordial power spectrum improves 
the total x 2 f° r the fit from 9803 to 9798. One expects that the effect of interest 
here would only make a small change (Ay 2 = 5) in the overall fit because it only 
affects a limited range of l values. Nevertheless, from the likelihood contours we can 
deduce mean value of A = 1.7 ± 1.5 with a maximum likelihood value of A = 1.5, 
and a mean value of fc* = 0.0011 ± 0.0004 h Mpc ” 1 
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k* 


Fig. 2. Constraints on parameters A and fc* from the MCMC analysis of the CMB power 
spectrum. Contours show 1 and 2 a limits. The horizontal axis is in units of [h Mpc -1 ). 


The values of A and fc* determined from from the CMB power spectrum relate to 
the inflaton coupling A and fermion mass m, for a given inflation model via Eqs. ([8]) 
and (|9|. 

A=\<j>4- 1 N\n*H- 1 ( 11 ) 


The coefficient A can be related directly to the coupling constant A using the 
approximatiorG3IHHIll for the particle production Bogoliubov coefficient 


|/3fc| 2 = exp 


—tt k 2 


IN A|</>* 


( 12 ) 


Then, 


This gives 


n* 


2 f°° 

— / dk p kp \/3k\ 2 
n Jo 


NX 3 / 2 

27r 3 


l ^| 3/2 


A = 



NX 5 / 2 1 

2V5ttV 2 


(13) 


(14) 

(15) 


where we have used the usual approximation for the primordial slow roll inflationary 
spectrum EH This means that regardless of the exact nature of the inflationary 
scenario, for any fixed inflationary spectrum 5ij(fc)U=o without the back reaction, 
we have the particle production giving a dip of the form Eq. with the parameter 
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A expressed in terms of the coupling constant through Eq. (151. Given that the 
CMB normalization requires Sn(k) | a=o ~ 10~ 5 , we then have 

A ~ 1.3NX 5/2 . (16) 

Hence, for the maximum likelihood value of A ~ 1.5, we have 

(1-1 -o.D 1 


N 


(17) 


A 5/2 A 5/2 ' 

So, A < 1 requires N > 1 as expected. 

The fermion particle mass m can then be deduced from m = N\(f>*. From 
Eq. (17) then we have m ~ </*/A 3 / 2 . For this purpose, however, one must adopt a 
specific form for the inflaton potential to determine <^>* appropriate to the scale fc*. 
Here we adopt a general monomial potential whereby: 




(18) 


for which there is a simple analytic relation 3 between the value of </>* and the number 
of e-folds A/”(&*) between when fc* exits the horizon and the end of inflation, i.e. 


Af(K) = 


77ipl" 


v(4>) 

V'W 


d<j) 


(19) 


implies 


/>* = VZaAfm p i . ( 20 ) 

For K = 0.0011 ±0.0004 h Mpc -1 , and k H = a 0 H 0 = (/i/3000) AIpc " 1 ~ 0.0002, 
we have AT—J\T* = /n(fc_y/fc*) < 1. Typically one expects A/’(fc*) ~ J\f ~ 50 although 
one can have the number of e-folds as low as J\f ~ 25 in the case of thermal inflation.® 
So, for a monomial potential with a = 2 we have 0* ~ 10 — 14 m p i. On the other 
hand, the limits on the tensor to scalar ration from the Planck analyst are more 
consistent with a = 1 (0* = 7—10 m p i), or even a = 2/3 (/>* =6 — 8 m p i). Hence, 
we have roughly the constraint, 

m ~ 6 - 14 ’ ( 21 ) 

well in excess of the Planck mass and independently of the number of degenerate 
species. 

As noted above, it is natural 10 to have such trans-Plankian massive particle in 
numerous extensions of the Standard Model such as supergravity and superstring 
theories. Such extra-dimensional theories generally contain a spectrum of particles 
with masses well in excess of the Planck mass. The extra-dimensions are compact 
and smaller than the three large spatial dimensions. It is, therefore, possible to 
dimensionally reduce the system to obtain an effective (3+1) dimensional theory 
that produces a tower of Kaluza-Klein (KK) statesE3H2l The mass of states in 
this tower is of order of the inverse size of the extra dimension. Since the extra 
dimensions are expected to have a size characteristic of the Planck length, these 
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KK states therefore have masses in excess of the Planck masses. Furthermore, these 
KK states can be nearly degenerate, with the level of degeneracy depending upon 
the geometrical structure of the compact space. It is, therefore, natural to deduce 
that a large number of nearly degenerate fermions existed during inflation with 
a mass well in excess of m p i, and that these particles couple to the inflaton field 
that drives inflation. For our purposes this degeneracy factor is described by the 
parameter TV, and in such theories N can easily of the order of 100. Hence, such a 
detectable signature in the CMB power is not unexpected. 


3. Matter Power Spectrum 

It is straight forward to determine the matter power spectrum. To convert the 
amplitude of the perturbation as each wave number k enters the horizon, 5n{k), 
to the present-day power spectrum, P(fc), which describes the amplitude of the 
fluctuation at a fixed time, one must make use of a transfer function, T{k) 21 which is 
easily computed using the CAMB codb^for various sets of cosmological parameters 
(e.g. f2, H 0 , A, f 1b)- An adequate approximate expression for the structure power 
spectrum is 

= - < 22) 

This expression is only valid in the linear regime, which in comoving wave number 
is up to approximately fc<0.2 h Mpc -1 and therefore adequate for our purposes. 
However, we also correct for the nonlinear evolution of the power spectrum. 22 

Figure [3] shows the matter power spectrum from the Wiggle-Z Dark Energy 
SurveTp ^ compared to the computed maximum likelihood power spectrum with 
and without the perturbation due to the resonant particle creation. Unfortunately, 
the perturbation is on a scale too large to be probed by the observed matter power 
spectrum. 


4. Low Velocity Cosmic Dark Flow and the CMB power spectrum 

An alternative way to explain the dip for £ = 10 — 30 has been pointed out by a 
number of authors 2 3 as due to an abrupt jump in the inflation-generating potential. 
In normal slow roll inflation, 


dV{(j>)/d(j) 
3 H 


(23) 


Hence, a large value for dV (c j> ) /d<f> can also produce a jump in the primordial power 
spectrim of Eq. [TJ In Gangopadhyay and MathewiPwe have explored this possibility 
with a potential of the form 


V(<p) = ^nup 2 + exp [(-¥>M)“] , 


(24) 
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Fig. 3. Comparison of the observed galaxy cluster function froirPU with the spectrum implied 
from the fits to the matter power spectrum with (solid line) and without (dashed line) resonant 
particle creation during inflation from Ref. [5] 


where, a is any odd number, but we adopt to a < 3 to ensure the renormalizibility. 
The parameter ipi then fixes the dip in the potential P For appropriate values of the 
potential parameters, a fit very similar to that of Figure [l] is attainable. 


4.1. Potential break and cosmic dark flow 

The particular relevance of our work, 8 however, is that a break in the inflation¬ 
generating potential not only explains the dip, but resolves another problem in 
early-universe cosmology, namely the magnitude of cosmic dark flow predicted in 
the M-theory landscape.^ 

The string landscape 24 is a complex space of a huge number of string vacua. 
These vacua differ from one another due to the complexity of RR fluxes and whether 
they are included or not and also the appearance of the gauge groups in the lower 
energy scale. A full analysis of the structure of the landscape is beyond the scope 
of present physics. However, there should be at least lO 500 to 10 6O ° possible vacua 
points in the landscape? 28 Starting from this, a model for the landscape was con¬ 
structed 2 ®^ as a lattice of vacua with a distribution of lattice energies. This land¬ 
scape lattice-like structure has two sectors SUSY and non-SUSY. Localization dis¬ 
connects these two sectors, while in the SUSY sector the landscape is viewed as a 
regular lattice, the non-SUSY sector exhibits an Anderson localization around each 
vacua site. Accepting that Supersymmetry is a broken symmetry, then the vacuum 
energy densities can vary widely being positive, negative or zero. 

All landscape vacua have equal probability to give birth to a universe. It has 
been shown in Refs. |25| - |27| that there are two kinds of baby universes: one kind will 
survive the back reaction due to quantum gravitational effects; and the other will be 
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terminal universes that will not survive the back reaction. Since all of the surviving 
universes were casually connected initially due to quantum unitarity, there will be 
an entanglement among all survivors. This entanglement should be imprinted in 
our survivor universe as well. 


4.2. Modified Friedmann Equation 

In the landscape scenario one must include the back-reaction due to the quantum 
entanglement with super horizon size wavelength modest 2 ' Then, the modified 
Friedmann equation becomes: 


H 2 = 


1 ,V(<p) 


V eff 


m [VM + 2m )F(b - v)]s m 


(25) 


where, 


F(b,V) 



m 2 Mp 

V 


)log{ 


b 2 M 2 
V ’ 


Ki¬ 


rn 2 b 2 Mp 

) exp (—3 y ) 


(26) 


Here, Mp is 8ttGn-, V{ip) is the inflation potential, b is the SUSY breaking scale 
and to 2 = V The interference length is given by: 


L 


2 

1 


a r ,m 2 . b m 2 H 1 

H l( 3Jf >'"Jf “ — ( P 



(27) 


4.3. ’Tilting’ gravitational potential 


The length scale L\ characterized the superhorizon size inhomogeneties. This non¬ 
local entanglement leads to a modification term (5<f> added to the gravitational po¬ 
tential. 


<F = $° + <5$ ~ <F°[1 


V(tp)F{b,V) r 
3 M 2 l T 1 


(28) 


where 4 >0 is the inflationary standard gravitational potential and <5$ is the ’tilt’. 
This tilt in the potential is responsible for the dark flowP3E 27 l 
For the quadruple anisotropy we geP 7 


AT 


■HI! 


Similarly, the effect on anisotropy for dipole is! 

AT, 4t t r r hw V(tp)F(b,V)^ 

T ldlp - 15 18Mp ’ 

The dark flow velocity (alid for L\ th) is then 2 ® 27 ! 


(29) 


(30) 


t AT M 

P — — cx( )\dip 


4?r / r h x , V(<p)F(b , V) 
15 L\ 18 M 2 


c 


(31) 
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Fig. 4. Dark velocity vs z plot with b ~ 10 10 Mp 


Where we have used 100 km s -1 for each 0.4 fiK anisotropy in converting the dipole 
temperature as in Ref. [27] 

We have revisecP the theory of Refs.EMIZ] f or a generic inflaton potential and 


applied it to the specific ’’break” potential described above in Eq. (24). Changing 


the inflaton potential to the new potential in Eq. (24) modifies the form of the 


correction term F(b,V). Since the limits on the SUSY breaking scale are quite 
stringent, 10 ~ 10 Mp < b < 10 ~ s Mp, the change in the modification mainly comes 
from the fact the inflaton potential is chaotic-like. This is an interesting case as 
most of the other possible potentials give much higher dark-flow velocity. 

In Figure [4] from Ref. 8] one can see the dark flow vs. redshift z for SUSY scale 
b ~ 10 ~ 10 Mp. This energy scale can be tested in the distant future in the proposed 
ILC. If we go to a higher breaking scale we get a higher value for dark flow. 

In Refs. [25]]27| a SUGRA like potential was employed of the form: 

<P 


V{<p) = U 0 exp(-A —) 


(32) 


Using this potential a ’dark’ flow of the order of 650-750 km s -1 was obtained. In 
Ref. [8] we have repeated the calculations of Refs. |25fj27] but with with the jump 
potential of Eq. (241. Adopting our model the resultant dark flow velocity diminishes 


considerably is illustrated in Figure [4] 


5. Observational Constraints on Dark Flow 

The above analysis of a lower landscape-generated dark flow velocity is of interest 
because there are now significant constraints on the bulk flow velocity. In Mathews 
et al! 2 ^ we have summarized recent constraints on the dark (or bulk) flow velocity, 
and completed a new analysis based upon an analysis of the SNIa redshift-distance 
relation. We made two analyses: one identifying the three dimensional Cartesian 
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velocity components; the other analyzing the cosine dependence on the sky of the 
deviation from Hubble flow. Fits were for z < 0.05 and z > 0.05 using both the 
Union2.1 and SDSS-II supernova surveys. We also studied® simulated data in which 
a bulk flow was imposed to determine whether the difficulty in detecting a bulk flow 
at high redsliift is due to uncertainty in the redshift-distance relation, confusion 
with peculiar velocities, or the absence of a bulk flow. We found 29 a bulk flow 
velocity of 270 ± 30 km s -1 in the direction ( l , b ) = (295 ± 30,10 ± 5)° in the 
Cartesian anaylsis, while the cosine analysis gave 325 ± 50 km s^ 1 in the direction 
(l,b) = (276 ± 15,37 ± 3)°, consistent with previous analyses. In the redshift bin 
z > 0.05, however, we found 29 only marginal evidence for a bulk flow velocity. We 
also found that the SDSS-II supernova data set has insufficient sky coverage to 
provide a meaningful result. 

The current observational situation is summarized in Ref. [29| as follows: it has 
been known since the 1980 ’d® that the local dipole flow extends well beyond the 
local super-cluster. This source was dubbed “the Great Attractor.” However, sub¬ 
sequent work in the 1990’s 31 has shown that the local flow extends at least to 
130 /i -1 Mpc. However, in Refs. 32 33 peculiar velocity fields were analyzed and 


found to be in agreement with a standard ACDM cosmology with no bulk flow, 
while in other work® an excess bulk flow was deduced. Moreover, there has little 
evidence of infall back toward the “Great Attractor” from material at larger dis¬ 
tances. Although there is recent evidencd® of a supercluster extending to a scale of 
~ 160 h -1 Mpc, there remains a need to search for evidence of a bulk flow at dis¬ 
tances well beyond ~ 160 h” 1 Mpc. Indeed, a recent analyst® 0 f the Cosmicflows-2 
galaxy peculiar velocity catalog has confirmed a bulk flow of 250 ±21 km s _1 out 
to a radius of 50 h -1 Mpc, and inferred a detectable bulk flow out to 200 h -1 Mpc. 

It has been proposed 37 ^® that detecting such dark flow on cosmic scales may be 
possible by means of the kinetic Sunyaev-Zeldovich (kSZ) effect. This is a distortion 
of the CMB spectrum along the line of sight to a distant galaxy cluster due to the 
motion of the cluster with respect to the frame of the CMB. Indeed, a detailed 
analysi^ 7 ®^ of the kSZ effect based upon the WMAP datcP seemed to indicate 
the existence of a large bulk flow velocity out to a distance of nearly 800 h” 1 Mpc. 
However, this is an exceedingly difficult analysis JUIEU and the existence of a dark 
flow has not been confirmed in a follow-up analyst 1 using the higher resolution 
data from the Planck Surveyor. The Planck results set a (95% confidence level) 
upper limit of 254 km s _1 for the bulk flow velocity and are consistent with no dark 
flow. It has, however, been argued® that the background averaging method in the 
Planck Collaboration analysis may have led to an obscuration of the effect. Indeed, 
a subsequent analysid® of the combined WMAP and Planck data has suggested 
evidence of a dark flow. Nevertheless, the Planck Collaboration upper limit is still 
consistent with as much as half of the observed CMB dipole corresponding to a 
cosmic dark flow. 

Independently of whether the dark flow has been detected via the kSZ effect, in 
view of the potential importance of this effect it is worthwhile to search for such 
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cosmic dark flow by other means. 

An analysis of peculiar velocity fields is an alternative means to look for the 
effect of a dark flow. In particular, an analysis of the distance-velocity relationship 
(and more importantly its residuals) well beyond the scale of the Local Group is 
needed to identify a cosmic bulk flow. A number of attempts along this line have 
been madeIMIE§MI 

For example, in Refs. 46 £7 the COMPOSITE sample of «4500 galaxies was 
utilized. For this set, distances were determined by the Tully-Fisher or “great plane” 
approach, plus a few by SN la distances. No dark flow was identified for distances 
greater than about 100 /i _1 Mpc. However, this combined data set might have large 
systematic uncertainties from the distant objects. Hence, a more carefully selected 
data set was subsequently analyzed.®* In that analysis the bulk flow magnitude was 
reduced to ~ 300 km s -1 and was only apparent out to about 80 h _1 Mpc. 

Galaxies in which a SN la has occurred provide, perhaps, the best alterna- 
tivd^2H55| because their distances are better determined. However there are fewer 
data available. There has been a wide variety in the attempts to find a bulk flow 
in SN la data sets. For example, Dai et al.J^ utilized the Union2 data set from the 
Supernova Cosmology ProjeclPi! and used a Markov Chain Monte Carlo (MCMC) 
search for the velocity and direction of the bulk flow. Colin et alP^did the same but 
used a different maximum likelihood method. Weyant et alP*^ utilized a unique data 
set as described in their paper. They then used a weighted least squares and a coef¬ 
ficient unbiased method to find the coefficient of the spherical harmonics describing 
the cosmic expansion. Feindt et alP*^ analyzed data from the Union2 compilation 
combined with that from the Nearby Supernova Factory!®’ They added a coherent 
motion into their cosmological model to search for a dark flow. 

The issue of sample sparseness has been considered in a number of works 
Weyant et alPS found that biases appear due to a non-uniform distribution or if there 
is power beyond the maximum multipole in the regression. Rathaus et al.^ also 
noted effects from large individual errors, poor sky coverage and the low redshift- 
volume-density, but they were still able to find a statistically significant bulk flow. 
Appleby et alP^ found a bias in galactic latitude and attributed it to the lack of 
data in the region of |6| < 20°. 

Our studjJ^Hl differed from that of the previous analyses in several key aspects. 
We made an independent analysis based upon two different approaches and two 
separate data sets. We began with a slight improvement of previous searchcfJMEEIl 
based upon the Union2 data set by analyzing ~ 600 galaxies with SN la redshifts 
from the Union2.1 supernova data setP^ We utilized a MCMC fit to the three 
Cartesian components of bulk flow velocity rather than the velocity magnitude in 
galactic coordinates. This approach had better stability near the Galactic pole. We 
then also analyzed the same data by searching for a deviation from Hubble flow with 
a cos 6 angular dependance on the sky. These studies established the robustness of 
these two complementary techniques for identifying the magnitude and direction 
of the bulk flow and confirmed previous detections of a bulk flow out to at least 
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z = 0.05. 

Having established the viability of the methods adopted,^ we then applied them 
for the first time to the large sample of (~ 1000) galactic redshifts and SN la 
distances from the SDSS-II survey HHOED However, we found that the analysis of the 
SDSS-ll data was severely limited by the paucity of data in the direction of the 
cosmic dipole moment. We establish^ that there is a detectable bulk flow at low 
redshifts, but at best a marginal detection for high redshifts. 

From simulated data sets, we deduced-® that the current uncertainty at high 
redshifts arises mostly from the current error in the distance modulus. We esti¬ 
mated® that with a sample like the Union2.1 data set, a detection would require 
both significant sky coverage of SNIa out to z = 0.3 and a distance modulus error 
reduction from 0.2 to <0.02 mag. However, a greatly expanded data set of 3xl0 4 
events might detect a bulk flow even with a typical distance modulus error of 0.2 
mag asmay be achievable with the next generation of large surveys like LSST. 
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